Critical fixed points in class D superconductors 
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We study in detail a critical line on the phase diagram of the Cho-Fisher network model separating 
three different phases: metallic and two distinct localized phases with different quantized thermal 
Hall conductances. This system describes non-interacting quasiparticles in disordered superconduc- 
tors that have neither time-reversal nor spin-rotational invariance. We find that in addition to a 
tricritical fixed point Wt on that critical line there exist an additional repulsive fixed point Wn 
(where the vortex disorder concentration Wn < Wt), which splits RG flow into opposite directions: 
toward a clean Ising model at W = and toward Wt- 
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The properties of quasiparticles in disordered super- 
conductors belonging to new symmetry classes^, in par- 
ticular, transitions between metallic, localized, or quan- 
tized Hall phasesS^^ , have been intensively studied. 
The symmetry class D may be realized in superconduc- 
tors with broken time-reversal and spin-rotation imparl- 
ances (as in d-wave superconductors with spin-orbit scat- 
tering). The associated changes in quasiparticle dynam- 
ics should be examined by energy transport, since nei- 
ther charge density nor spin are conserved. In this Brief 
Report we present a detailed study of the crtitcal line 
on the phase diagram for a model first introduced by 
Cho and Fisher (CF)£ (see a detailed description be- 
low) , which has a particularly rich phase diagram in two 
dimensions 5 ^. 

We employ advanced numerical calculations, proposed 
in Ref. [fi and described in detail in&2 to overcome 
round-off errors in calculations of renormalizcd localiza- 
tion lengths. We also apply an optimization algorithm 
to determine both critical exponent and critical energy. 
These results allow us to determine the tricritical point 
where three phases meet and to study dependence of crit- 
ical exponent for the insulator-to-insulator transition on 
the width of the system. Our results suggest the exis- 
tence of two fixed points. 

The original network models was proposed to describe 
transitions between plateaux in the quantum Hall effect 
(QHE). In the model flux probabilities move along unidi- 
rectional links forming closed loops in analogy with semi- 
classical motion of electrons on contours of constant po- 
tential. Scattering between links is allowed at nodes in 
order to map tunneling through saddle point potentials. 
Propagation along links is described by diagonal matri- 
ces with elements in the form exp(i</>). Transfer matrix 
for one node relates a pair of incoming and outgoing am- 
plitudes on the left to a corresponding pair on the right; 
it has the form 
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where e is a dimensionless relative distance between the 
electron energy and the barrier height. It is easy to 
see that the most "quantum" case (equal probabilities 
to scatter to the left and to the right) is at e = 0. 



Numerical simulations on the network model are per- 
formed on a system with fixed width M and periodic 
boundary conditions in the transverse direction. By mul- 
tiplying transfer matrices for N slices and then diagonal- 
izing the resulting total transfer matrix, it is possible to 
extract the smallest Lyapunov exponent A (the eigenval- 
ues of the transfer matrix are exp(AiV)). The localization 
length £a/ is proportional to 1/A. Renormalized local- 
ization lengths for different system widths and different 
energies satisfy a one-parameter scaling 
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which provides dependence of the thermodynamic local- 
ization length £ on energy e. 

For a D class symmetry a Bogoliubov - de Gennes 
Hamiltonian is written in terms of a Hermitian matrbi 
The corresponding time evolution operator is real, and, 
the generalized phase factors are, therefore, O(N) ma- 
trices for a model in which iV-component fermions prop- 
agate on links. We study the case N = 1 with phase 
factors ±1. There are three models: random bond Ising 
model— supporting two different localized phases, un- 
corrected 0(1) modeUS, where phases on the links are in- 
dependent random variables and all states are extended 5 , 
and the model first introduced by Cho and Fisher (CF) 6 
where scattering phases with the value tt appear in cor- 
related pairs. Each model has two parameters: the first 
one is a disorder concentration W, such that there is 
a probability W to have a phase 0, and a probability 
(1 — W) to have a phase it on a given link. The sec- 
ond parameter is an energy e describing scattering at the 
nodes. For the CF model, the phase diagram (updated 
version of which is presented in Fig. 1) in the e-W plane 
has three distinctive phases: metallic, and two insulating 
phases characterized by different Hall conductances. The 
sensitivity to the disorder is a distinctive feature of class 
D. 

In the CF model the disorder is introduced only at the 
nodes, allowing for the offdiagonal elements of Eq. (1) to 
be multiplied by ±1 (disorder probablity W is the proba- 
bility of that factor to be —1). In our previous worki^ we 
have studied a CF phase diagram and its critical expo- 
nents far from e = 0. On the other hand, the critical line 
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itself at e = is of the particular interest to us. One of 
the possibilities discussed in Ref. [H, where a quasipar- 
ticle density of states for the CF model was studied, is 
the existence of the second fixed point (a repulsive one) 
on the critical line separating two insulating phases at 
the disorder smaller than the one at the tricitical point, 
where the critical line splits. It has been suggested 16 that 
we can address this question using our optimization pro- 
cedure studying the analogue of RG flow for the critical 
exponent. Indeed, if one belives the scenario of two fixed 
points: the repulsive one at Wn (the subscript N sug- 
gested in Ref. [TH to underline the similarity to the role 
of the Nishimori point in the RBIM) and the tricritical 
one at Wt > Wn, then the one-parameter scaling for the 
disorder W < Wn should produce the critical exponents 
which tend to the value v = 1 of the pure Ising model as 
one uses only large system widths. In contradistinction, 
for W > Wn the repulsive fixed point should push the 
crtitical exponent to flow toward the metallic fixed point. 

We have performed numerical calculations at small 
fixed values of W for the system widths M = 
16,32,64, 128 and various energies e. In those cases an 
obvious critical energy is e = (the most " quantum" case 
explained above) . We have been determining the critical 
exponents v (£ ~ e~ v ) to fit all the data onto one curve 
by applying a special optimization program which checks 
different critical exponents and chooses the optimal one. 
We have carried out the analysis, first, for all data points 
for all system widths, second, without M = 16 data, 
third, without M = 16, 32 data. It turns out that, in- 
deed, this procedure shows that the critical exponent 
flows to the predicted pure Ising model value v = 1 as we 
omit small system widths data (the analogue of RG flow) . 
As a typical example we present here the flow of the crit- 
ical exponent v for the disorder W = 0.04: v = 1.34 
for M = 16,32,64,128, v = 1.11 for M = 32,64,128, 
v = 0.97 for M — 64, 128. The largest value of disorder 
for which this flow persists is 0.1. To demonstrate the 
flow in a conclusive way at this disorder value, the data 
for the largest system width M = 256 was necessary: the 
critical exponent flows as 1.7 — ► 1.6 — > 1.2 — » 1.1. We, 
therefore, identify Wn = 0.1. In a full agreement with 
the scenario suggested^ for W > Wn, the critcal expo- 
nent flows in the opposite direction, e.g. for W — 0.12 
the critical exponent is almost constant « 1.7, indepen- 



dent on the system widths used. In order to be sure that 
we are still on the critical line we allow the oprimiza- 
tion program to look for the critical energy as well. It 
is very conclusive that up to the disorder value 0.14 the 
critical energy is about 10~ 5 which supports the single 
critical line. For W > 0.14 the optimization program 
immediately produces small but finite values of the crit- 
ical energy (as shown on Fig. 1). We thus identify the 
triritical fixed point Wt = 0.14. 

To summarize, we have studied in detail the critical 
line on the phase diagram of the CF model. We have 
used the RG flow of the critical exponent and determined 
two fixed points in agreement with one of the scenarios 
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FIG. 1: Updated phase diagram for a CF model with metallic, 
insulating and quantized Hall phases. 



suggested in Ref. Qjl: the repulsive fixed point at disorder 
Wn = 0.1 and the tricritical fixed point Wt = 0.14. 
Obviously, we cannot rule out completely that for larger 
system widths the RG flow can reverse as was founcU^. 
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